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ABSTRACT 

Direct TV-body calculations are presented of the early evolution of exposed clusters 
to quantify the influence of gas expulsion on the time-varying surface brightness. By 
assuming that the embedded OB stars drive out most of the gas after a given time 
delay, the change of the surface brightness of expanding star clusters is studied. The 
influence of stellar dynamics and stellar evolution is discussed. 

The growth of the core radii of such models shows a remarkable core re-virialisation. 
The decrease of the surface mass density during gas expulsion is large and is only 
truncated by this re-virialisation process. However, the surface brightness within a 
certain radius does not increase noticeably. Thus, an embedded star cluster cannot 
reappear in observational surveys after re-virialisation. This finding has a bearing on 
the observed infant mortality fraction. 

Key words: methods: TV-body simulations - galaxy: open clusters and associations: 
general 



1 INTRODUCTION 

Since stars form in dense clumps within molecular clouds, 
very young clusters are observed to be embedded in gas. 
Typically 1/3 of the initial cluster-forming gas mass is con- 
verted into stars (Lada & Lada 2003; Machida & Matsumoto 
2011). Following star formation, and especially after forma- 
tion of OB stars, most of the residual gas is driven out, so 
that there is nearly no gas observed in populations older 
than a few million years. While the gas ejection process 
begins, self-gravitation decreases, resulting in a rapid ex- 
pansion of the young star cluster (Goodwin 1997, Good- 
win & Bastian 2006, Boily & Kroupa 2003, Baumgardt & 
Kroupa 2007). In this contribution we investigate cluster 
formation and evolution during this process, focusing on 
the surface brightness of star clusters responding to residual 
star-forming gas expulsion. We want to answer the question 
of how the surface brightness is influenced by a) stellar dy- 
namics and b) stellar evolution, and how both relate to each 
other. 

In the context of gas expulsion, numerical A-body com- 
putations show a re-virialisation process taking place after 
a few million years in the centre of the system. While the 
cluster outskirts continue to expand, the centre contracts 
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(Kroupa, Aarseth & Hurley 2001), leading to an increase of 
the density in the central region. Since at any given age, the 
surface mass density has an impact on the surface bright- 
ness, this process is of special interest. We ask what role 
this effect plays in the evolution of the surface brightness 
and whether its impact is sufficient to let the cluster core 
disappear under a given detection limit owing to its expan- 
sion, then re-appear after re-virialisation. 

We have performed direct A-body computations, in- 
cluding stellar evolution, with initial conditions similar to 
Kroupa et al. (2001), with the emphasis being put on the 
time-evolution of the cluster surface brightness and sur- 
face mass density. The residual gas is modelled by a time- 
dependent background potential. 



2 THE CODE 

The main aspect here is the impact of gas expulsion on 
the early evolution of star clusters. We choose Nbody 6 
(Aarseth 2000) to perform these direct A-body simulations 
since this code already offers features that enable us to ac- 
count for gas expulsion. In doing so, a time-varying analyt- 
ical background potential is added which represents the gas 
and therefore its gravitational interaction with the stellar 
population. We have to make this assumption to simplify 
the complex hydrodynamical processes in the gas compo- 
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nent. However, it has been concluded by Geyer & Burkert 
(2001) that this approximation is physically realistic. 

In addition to the gas potential, a standard solar- 
neighbourhood tidal field is adopted. Also, stellar evolution 
is taken into account, because it is important when deter- 
mining the luminosity of stars. We did not include pre-main 
sequence evolution as it is not part of the Nbody 6 pack- 
age yet, and incorporating it into Nbody 6 is a very major 
effort. This means that the star clusters would initially be 
brighter than we compute, but also that they would fade 
even more strongly. The brightening of some of the clusters 
listed in Tab. 2 are therefore upper limits, since pre-main se- 
quence evolution would reduce this brightening because the 
stars fade towards the main sequence. The brightening effect 
due to revirialisation is therefore even less observable than 
estimated by our present models such that clusters cannot 
reappear after revirialisation. 

2.1 Analytical gas model 

We assume that both the stars and gas obey a Plummer den- 
sity profile (Plummer 1911) with equal half-mass radii. The 
star formation efficiency (SFE) at the onset of gas expulsion 
is SFE = 1/3, that is, the gas mass is twice the initial stellar 
mass M Bt (t < to), 



M e (t < to) =2M st (t < to) . 



(1) 



A star with a position r relative to the centre of the 
cluster experiences an acceleration from the background po- 
tential, 



GM s (t) 
(r»+r; iig )3/ 



(2) 



where M g (t) is the time- varying mass and r p i jg is the Plum- 
mer radius of the cluster-forming region. G is the gravita- 
tional constant. 

Gas expulsion is modelled by decreasing the gas mass, 
the Plummer radius r p i !g remaining unchanged. The gas 
mass M g (t) - and accordingly the background potential - 
is assumed to stay constant until a time-delay to = 0.6 Myr. 
After this time delay it decreases following 



M s (t) = 



M e {t < to) 



(3) 



l + (t-te)/r' 

where r is the gas expulsion time scale, which depends on 
the velocity of the heated ionized gas (lOkms -1 , e.g. Hills 
1980), i.e. 

r pi.g 



lOpc/Myr 



(4) 



We stress that the impact of the gas expulsion time scale 
depends on r expressed in units of the cluster-forming region 
crossing time, t CT : 



8r 3 . 

vir 

GM st 



where the virial radius obeys 



- 1£ 
r vir - ^r p i 



(5) 



(6) 



in a Plummer model with a Plummer radius r p \. Gas re- 
moval can be expected to be very rapid, and faster than 



the dynamical crossing time of the cluster-forming region 
(Whitworth 1979). In the present calculations, r takes the 
following values 



t £ {0.001, 0.15, 0.3, 0.6, 1.2} t c 



(7) 



Table 1 gives t CT in Myr. The time-delay to is adopted from 
Kroupa et al. (2001) and corresponds to the gas-confinement 
time until the HII region erupts. At a time t = to + t, 
half of the initial gas mass is driven out. In that respect, r 
differs from the gas expulsion time scale adopted in Geyer 
& Burkert (2001) and Parmentier et al. (2008), where it is 
defined as the time when gas removal is complete. 



Alternative analytical model 

One alternative model representing the gas would be a 
Plummer model with constant mass and time-dependent 
Plummer radius r p \, s (t), e.g. 



r (o) 
'pi,g' 



1 + 



- 1 



1/2 



* < to 
t ^ to , 



(8) 



where r^j is the Plummer radius at t < to- This expres- 
sion is derived from the condition to have, at the time t, 
the same gas mass left within a sphere about the centre 
of the system with radius = r p i jg (0) as in the first 

gas model. In Fig. 1, the corresponding density and poten- 
tial are plotted for t £ {to, to + t, to + 10r} (dashed line). 
Compared to the model based on Eq. 2 we use in our cal- 
culations (solid line), the gas mass is not just reduced but 
the gas is driven outwards. As a consequence, in the model 
with variable Plummer radius the gas is removed faster from 
the central region (see Fig. 1, first panel), and the potential 
remains deeper (second panel) so that cluster will disrupt 
more slowly as long as the Plummer radius of the gas is 
smaller than the tidal radius. The maximum of the reject- 
ing force, —d(f>(r)/dr, resulting from the gas potential can 
be found at 



c(*) 



r pl,g(*) 

V2 



(9) 



Thus, increasing the Plummer radius, r p i :g (t), means in- 
creasing the radius r max (i) at which the rejecting force is 
maximal, resulting in a larger, less dense core. 

In the present work, we use the decreasing mass model, 
because we assume the star cluster models to be initially not 
mass-segregated. In this case, it is more natural to remove 
the gas uniformly, following the distribution of the OB stars. 



3 INITIAL CONDITIONS 

Five models with different initial properties (total mass and 
half-mass radius) are calculated. Model III is based on model 
A of Kroupa et al. (2001). The properties of all models are 
listed in Tab. 1. 

As for the tidal field, a Milky Way potential is adopted. 
The cluster has a distance to the Galactic centre of 8.5 kpc 
and moves with a velocity of 220 km s -1 on a circular orbit 
through the Galactic disc. 
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Figure 1. Comparison of two different analytical gas expulsion 
models. Both models assume that the initial gas distribution (at 
t < to) follows a Plummcr model. The solid line corresponds to a 
model where the gas mass M g (i) is decreased with time (Eq. 3). 
In the alternative model (dotted line), the total gas mass is as- 
sumed to be constant, and the Plummer radius r p i ig is increased 
with time after t > to- The first panel shows the assumed gas 
density p g (r), while the second panel presents the resulting po- 
tential g (r), and the third panel shows the resulting radial force, 
mr = — d</> g (r)/dr. The plots compare the gas models at different 
times, to = to , ti = to ~f~ T, and tio = to + lOr. 



3.1 The initial stellar population 

The stars are initially distributed in mass according to the 
canonical two-part power-law IMF (Kroupa 2001). The dis- 
tribution function £(m) is given by 

£{m)cxm- a , (10) 



Table 

here is 


1. Initial properties of the models, 
model A of Kroupa ct al. (2001). 


Note that model III 




M st [Mq] 


N 


(m) [M ] 


Ro.s [pc] 


tcr [Myr] 


I 


787 


1343 


0.59 


0.263 


0.15 


11 


787 


1357 


0.58 


0.45 


0.33 


III 


3700 


6382 


0.58 


0.45 


0.15 


IV 


18893 


32578 


0.58 


0.45 


0.067 


V 


18893 


32122 


0.59 


0.77 


0.15 



This table gives an overview of the models we use. M B t is the total 
stellar mass before gas expulsion onset, TV the number of stars, 
(m) the mean stellar mass, Ro.5 the half-mass radius, and t cr the 
crossing time. The gas expulsion time scale r (see Section 2.1) 
takes values of 0.001t cr , 0.15t cr , 0.29 t cr , 0.59 t cr , and 1.18t cr in 
each model. 

where 

f 1.3, 0.08 s: m/M < 0.5, . . 

Q ~ \ 2.3, 0.5 s: m/M < 100. ( ' 

£(m)dm is the number of stars with mass m £ [m, m + dm]. 
For computational feasibility, there are no primordial bina- 
ries. This is a valid approximation, because they would not 
have an important effect on the surface brightness. 

3.2 Cluster model 

Like the gas, the density distribution of the stellar popula- 
tion at t < to with stellar mass M st is assumed to follow a 
Plummer model with a half mass radius, -Ro.5, that is iden- 
tical to that of the gas. The properties are listed in Tab. 1. 
The models are initially not mass-segregated. 



4 ANALYSIS STEPS 

For every time step, the Lagrange radii as well as the surface 
brightness are determined using the mass, position and to- 
tal luminosity of the stellar ensemble. The density centre is 
calculated by the method of Casertano & Hut (1985). Mass 
loss due to stellar evolution is already taken into account by 
the simulation. 

In order to analyse the time-evolution of the surface 
brightness, the average surface brightness within a constant 
radius (1, 3, 5, 10, and 15 Myr) is plotted vs. time. 

In order to compute the surface brightness of a finite en- 
semble of stars the stellar luminosity evolution as a function 
of the initial stellar mass is required. We use the spectral 
evolution code Pegase by Fioc & Rocca-Volmerange (1997, 
1999). This code allows a user-defined multi-part power law 
IMF. In order to extract the evolution for a single star with 
mass mo, we define a Dirac-IMF, £(m) = S(m — mo), and 
approximate it with a single-part power law, £(m) oc mT a , 
with slope a = between (1 — 0.001)mo and (1 + 0.001)mo 
and zero otherwise. The masses mo are the same as the 
masses for which stellar evolution models are implemented 
into Pegase. A two-dimensional table (the stellar luminos- 
ity as a function of the initial stellar mass and the age of the 
star) for the passbands FUV, B, U, and V are constructed. 
A luminosity for an arbitrary initial stellar mass and age 
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Figure 4. This plot illustrates the effects of stellar evolution and 
stellar dynamics on the cluster surface brightness. It shows the 
total FUV luminosity of all stars located within the innermost 
1 pc radius (dashed line, in ergs" ) and the FUV surface 

brightness of a projected area with the same radius (solid line, 
in ergs - 1 A J pc~ 2 ) on the left jy-axis. In order to relate the sur- 
face brightness evolution to the re-virialisation process, the 20% 
Lagrange radius (LR, grey line) is shown on the right y-axis. 

is obtained by two-dimensional linear interpolation of the 
luminosity table. 



5 RESULTS 

As we analyze the results of the performed computations 
qualitatively, only model III is presented in detail, because 
the remaining models show qualitatively the same results. 

The computational results are shown in Figs. 2 and 3. 
As a check of consistency, the evolution of the Lagrange 
radii is compared between model A of Kroupa et al. (2001) 
and model III here. They show the same characteristics: 
the expansion during and immediately after gas expulsion, 
followed by re-virialisation. 

In model III, the Lagrange radii, the surface mass 
density, and the surface brightness in the four passbands 
specified above are plotted vs. time (six columns in to- 
tal). The first column shows the Lagrange radii which doc- 
ument the expansion history of the cluster. Columns 2-6 
contain the time-evolution of the average surface mass den- 
sity/brightness within constant radii (1, 3, 5, 10, and 15 pc). 
Each row presents a computation of model III with a differ- 
ent gas expulsion time scale r. 

Figure 4 shows the 20% Langrage radius, the total 
cluster luminosity and the average FUV surface brightness 
within a projected 1 pc radius circle around the density cen- 
tre of the system. On the basis of this figure, we want to 
discuss the interaction of cluster expansion and total stellar 
luminosity evolution and, accordingly, the dependencies of 
the surface brightness on stellar dynamics and stellar evolu- 
tion. 

While the cluster starts to expand after the time to 
has passed, the total luminosity of all stars increases clearly 
until t « 2 Myr. It reaches its maximum near the core re- 
virialisation phase, depending on the model. From then on, 
the total stellar luminosity (dashed line) decreases distinctly 



Table 2. Re-increase of the surface brightness due to core re- 
virialisation in the V-band. 



Model 


T [tcr] 


SB [mag] 


SB min [mag] 


SB max [mag] 


III 


0.6 


-8.5 


-7.41 


-7.81 


IV 


0.001 


-9.1 


-8.3 


-8.7 


IV 


1.2 


-9.1 


-8.5 


-8.8 


V 


0.001 


-9.5 


-6.8 


-8.0 


V 


0.3 


-9.5 


-8.4 


-8.7 


V 


0.6 


-9.5 


-7.7 


-8.7 


V 


1.2 


-9.5 


-7.8 


-8.5 



The table gives an overview of all calculations with a positive re- 
increase of the surface brightness within the inner 1 pc due to core 
re-virialisation. Only 7 of 25 calculations show this behaviour. 
Here, t is the gas expulsion time, SBo the initial surface brightness 
at t < trj , SB m i n corresponds to the minimum surface brightness 
before re-virialisation at t m i n , and SB max to the maximum surface 
brightness at t m ax after re-virialisation. 

as massive stars evolve. Although the contraction of the core 
can be weakly detected as a bump in the surface mass den- 
sity (see Fig. 2, middle column), it is covered in the surface 
brightness by luminosity changes due to stellar evolution so 
that stellar dynamics play only a minor role in that context. 
Therefore, we can conclude that the re-virialisation in the 
core does not result in a significant re-increase of the central 
cluster surface brightness. 

To quantify the latter finding in all models (five models, 
each with five different values of r, i.e. 25 calculations in to- 
tal), we do the following analyses. First, look at the surface 
mass density and detect the re-increase of the surface mass 
density within the inner 1 pc which results from the core 
re-virialisation. If there is a noticeable increase, remember 
the times at which the surface mass density reaches the cor- 
responding local minimum, t m i n , and maximum, t max , and 
determine the surface brightness changes (in the V-band) 
between £ m i n and i m ax- Table 2 shows the resulting surface 
brightness changes of all calculations in which such a re- 
increase is detectable. 7 of 25 models match the require- 
ments. In model I and II, not a single calculation matching 
these criteria was found. In model III, r = 0.6 t CI results in a 
weak re-increase of the surface brightness. Model IV has also 
two calculations with a weak increase and model V shows 
a visible effect in four calculations. However, compared to 
the later surface brightness decrease due to stellar evolution, 
the detected changes, SB max — SB m i n , are fairly small and 
it seems unlikely that they lead to the addressed effect of 
falling below and (re-)exceeding the detection limit. 

Apart from our main question we notice some general 
dynamical aspects of gas expulsion when comparing sim- 
ulations with the same initial conditions and different gas 
expulsion rates. The expansion of the 90% Lagrange ra- 
dius docs not depend strongly on the gas expulsion time 
scale, although this parameter covers a range of three or- 
ders of magnitude. In contrast, the expansion of the inner 
Lagrange radii is very sensitive to the gas expulsion time 
scale. The quicker the gas expulsion, the faster the expan- 
sion of the central part of the cluster and accordingly of the 
inner Lagrange radii (first column in Fig. 2). This is due to 
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Figure 2. Model III These plots show the main results of the simulations. The left column presents the Lagrange radii (within each 
panel bottom-up: 2%, 5%, 10%, 20%, 30%, 40%, 50%, 80%). From bottom to top, the gas expulsion rate slows: the smaller r, the quicker 
gas expulsion. The solid thick curve describes the mass of the residual gas, which is represented by the second y-axis and is given in units 
of the total stellar mass before gas expulsion onset. In the middle and right columns, the evolution of the surface mass density and the 
surface brightness within the inner 1, 3, 5, 10 and 15 pc (solid curves, top to bottom) are presented. 
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Figure 3. Model III Complcmental to Fig. 2, this plot shows the surface brightness in the U, V, and B bands. 



the effective crossing time, which gets longer with increasing 
radius and is shortest in the cluster core. Consequently, in 
the case of slow expulsion re-virialisation influences the sur- 
face mass density less strongly, because the cluster remains 
bound anyway for a long period. For short expulsion time 
scales on the other hand, the core expands fast and the core 
has only little chance to re-virialise. For the assumed SFE, 



one can observe re-virialisation clearly only for a range of 
gas expulsion times r, depending on the model. 



The surface mass density plots show statistical fluctua- 
tions in the central region due to the small number of stars 
(the noisy shape of the upper line of the second row of Fig. 2 
after a few million years). The noisy shape in the surface 
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brightness plots on the other hand is mainly caused by the 
evolution of individual stars. 



6 CONCLUSIONS 

In the present work, we investigate the evolution of the sur- 
face brightness of star clusters under the assumption that 
the residual gas of the cluster-forming region is driven out 
by evolving OB stars, leading to a strong expansion of the 
total cluster. Five different models of star clusters with dif- 
ferent sizes and masses are numerically evolved, every model 
with a SFE of 1 /3 and with five different gas expulsion time 
scales r. Thereby we compare the influences of stellar dy- 
namics and stellar evolution on the surface brightness. Based 
on the expansion histories the question came up whether the 
impact of the re-virialisation process is sufficient to let the 
cluster core disappear under a given detection limit owing to 
its expansion, then re-appear after re-virialisation. From our 
computations we conclude that the surface brightness is not 
changed significantly by this process, because it is clearly 
dominated by stellar evolution. 

Independent by this finding, we can draw conclusions 
of the expansion behaviour. We find that the expansion of 
the inner parts of a star cluster is more sensitive to the gas 
expulsion time scale r than the outer parts are. While the 
outer Lagrange radii are only weakly responsive to changing 
t, the inner Lagrange radii react very sensitively. E.g. the 
core expands moderately in the case of slow expulsion, while 
a quick expulsion on the other hand can lead to a rapid 
disruption of the core. 
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